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| We have studied mixed states in the system of three qubits with the property that all their 

■ partial transposes are positive, these are called PPT states. We classify a PPT state by the ranks 
of the state itself and its three single partial transposes. In random numerical searches we find 
entangled PPT states with a large variety of rank combinations. For ranks equal to five or higher 
we find both extremal and nonextremal PPT states of nearly every rank combination, with the 
restriction that the square sum of the four ranks of an extremal PPT state can be at most 193. We 
have studied especially the rank four entangled PPT states, which are found to have rank four 
for every partial transpose. These states are all extremal, because of the previously known re- 
sult that every PPT state of rank three or less is separable. We find two distinct classes of rank 
4444 entangled PPT states, identified by a real valued quadratic expression invariant under local 
SL(2, (D) transformations, mathematically equivalent to Lorentz transformations. This quadratic 
Lorentz invariant is nonzero for one class of states (type I in our terminology) and zero for the 
other class (type II). The previously known states based on unextendible product bases is a non- 

. generic subclass of the type I states. We present analytical constructions of states of both types, 

general enough to reproduce all the rank 4444 PPT states we have found numerically. We can not 

■ exclude the possibility that there exist nongeneric rank four PPT states that we do not find in our 
| random numerical searches. 

t-H . 

(N| ■ PACS numbers: 03.67.Mn, 02.40.Ft, 03.65. Ud 

1 Introduction 

As one of the most fascinating features of quantum mechanics entanglement has been intensively 
studied in the last decades 0]. The Bell inequalities [2] turned the philosophical discussion between 
Einstein and Bohr into a subject for experimental investigations (3l 01 \5\ HI . The Bell inequalities 
are deduced from the hypothesis of local realism, and apply to statistical correlations in composite 
systems with two subsystems. One basic weakness of all such experiments is the inherent statistical 
uncertainties of the observed correlations. 

The striking new feature of the experiments proposed by Greenberger, Home, and Zeilinger ||7] 
and Mermin [8], in composite systems with three or more subsystems, is that the correlations to be 
tested are absolute and no longer statistical. One single observation is sufficient to demonstrate that 
quantum mechanics violates local realism. Experiments of this kind have also been made (9). 

A central problem in the study of entanglement is how to determine, theoretically or experimen- 
tally, whether a state in a composite system made up of two or more components is entangled or 



1 



separable. The answer is simple for a pure state: it is entangled if it is not a tensor product. The 
separability problem for mixed states, on the other hand, has proven to be highly nonuivial iTTOl and 
does not yet have a solution which is satisfactory for practical use. The simplest and best known 
condition for the separability of a mixed state is the Peres criterion ITTTI . which states that a separable 
state remains positive semidefinite under partial transpositions. This necessary condition is in general 
not sufficient, but applying it requires very little computational effort and the separability problem 
is therefore in essence reduced to determining whether mixed states with positive partial transposes 
(PPT states) are entangled or separable. 

A partial transposition is a mathematical operation that Uansposes one or more factors in a tensor 
product, for example, 

(pi ® P2® Pzf 2 = Pi® P2 ® P3 • (1) 

The point is that it is a well defined operation even for a matrix which is not a tensor product (see 
Appendix [A]). A Hermitian matrix representing a mixed state of an n-partite system may appear in 
2 n different versions related by partial transpositions, but these can be separated into two sets where 
every matrix in one set is the total transpose of a matrix from the other set. Since the total transpose 
T preserves eigenvalues, half of the 2™ partial transposes are superfluous when we evaluate whether 
or not a mixed state p is a PPT state. 

In the present article we consider PPT states in a system of three qubits. In this system it is 
sufficient to watch the eigenvalues of p, p Tl , p T2 , and p T ' A , where Tj is the partial transposition with 
respect to subsystem i. It is necessary to require that all these four are positive, since it is quite possible 
for three of them to be positive, while the fourth one is not. The four remaining partial transposes are 
obtained from these by a total transposition. For example, T1T2 = T3T, because T = T1T2T3, partial 
transpositions commute, and T { 2 is the identity operation. 

We write V for the set of all unnormalized mixed states, and T>\ for the set of all mixed states 
normalized to unit trace. Similarly, we write S or S\ for the set of all separable states, and V or V\ 
for the set of all PPT states, unnormalized or normalized. All of these are convex sets. 

An extremal point in a convex set is a point which is not a convex combination of other points in 
the set. The sets of normalized states, T>\ , S\ , and V\ are compact and hence completely described by 
their extremal points, in the sense that all nonextremal points may be written as convex combinations 
of the extremal points. The extremal points of T>\ are the pure states, and the extremal points of S\ 
are the pure product states. The inclusions 

5iCPiCDi, (2) 

together with the fact that the pure product states are extremal points in both S\ and T>\, imply that the 
pure product states are also extremal points of V\ . It is easy to prove that all pure nonproduct states 
are not in V\ and therefore entangled. The pure product states are the only PPT states of rank one, 
since only pure states have rank one. Because S\ and V\ are not identical in the three qubit system, V\ 
must have extremal points of rank higher than one giving rise to all the entangled PPT states. These 
extremal entangled PPT states are almost completely unknown, and that situation motivates our study 
presented here. 

Previous studies of multipartite entanglement have been mostly concerned with pure states. Diir, 
Vidal, and Cirac lfl2l classified pure states in the three qubit system into six equivalence classes, based 
on the type of entanglement possessed by a state. A pure state is of the form ipi/j' with ip G C 8 . Two 
(unnormalized) vectors ip and <fi are considered equivalent if 

(j) = (Fi ® V 2 <g> F 3 ) %l> (3) 
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with Vi G SL(2, C). One class contains the separable (unentangled) pure states where ?/>isa2x2x2 
dimensional product vector. Three other classes contain the biseparable pure states where ip is a 
product vector in one of the three splittings into one system of one qubit and one system of two qubits. 
All states in these three classes have only bipartite entanglement. The last two classes contain states 
with two inequivalent types of genuine tripartite entanglement. There is the W class, exemplified by 
the unnormalized state, in Dirac notation, 



W) = |100) + |010) + |001) 



(4) 



Finally there is the GHZ class, exemplified by the Greenberger-Horne-Zeilinger state 

| GHZ) = |000) + | 111) . 



(5) 



In matrix notation, with 



the W and GHZ states are 
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The corresponding classification scheme for mixed states in the three qubit system was introduced 
by Acm et al. Ifl3l . The scheme involves four convex sets with an onion structure where each set is 
defined by including increasingly larger sets of pure states. The innermost set is the set S of all 
separable states, consisting of all states that are convex combinations of pure product states. The 
second set B includes all the biseparable pure states. The third set W includes also pure states with 
W entanglement. The fourth set includes the last remaining class of pure states, the GHZ entangled 
states. Because it includes all pure states the fourth set is V, the set of all mixed states. The authors 
conjecture that all entangled PPT states are members of the third set W. 

Karnas and Lewenstein lfl4l proved that all states of ranks two and three in the three qubit system 
are separable. They applied the range criterion of entanglement for the three qubit system as a special 
case of the 2 x 2 x N dimensional system. Bennett et al. Ifl5l [161 introduced the PPT states based 
on unextendible product bases (UPBs) as examples of entangled PPT states in the three qubit system. 
These three qubit UPB states have rank four. They are entangled because there is no product vector in 
the range of a UPB state, and they are extremal PPT states because no entangled states of lower rank 
exist. 



Sorting states into equivalence classes 

We want to classify the three qubit states into SL <g) SL <g> SL equivalence classes, defining two unnor- 
malized density matrices p and a to be equivalent if 

a = (Vi <g> v 2 ® v 3 ) p (yi ®y 2 ® y 3 ) t , (8) 
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with Vi € SL(2, C). This definition is useful because equivalent density matrices have the same 
entanglement properties, although quantitative measures of entanglement will be different. Qualita- 
tive properties will be the same, such as tensor product structure of pure states, decomposition of 
mixed states as convex combinations of pure states, rank and positivity of states and all their partial 
transposes, and so on. 

The relation between the group SL(2,C) and the group of continuous Lorentz transformations is 
well known, and is reviewed here in Appendix |B] From a density matrix in the three qubit system 
we define one quadratic and four quartic real Lorentz invariants, so called because they are invariant 
under SL (g> SL (g> SL transformations as in equation ([8]). They are also invariant under partial transposi- 
tions, because a partial transposition may be interpreted as a parity transformation, which is a discrete 
Lorentz transformation. 

This means, for example, that the ratio between one quartic Lorentz invariant and the square of 
the quadratic invariant will have the same value for all the states in one equivalence class and all their 
partial transposes. Taking the ratio between Lorentz invariants is necessary in order to cancel out any 
normalization factor in the density matrix. If two density matrices are not in the same equivalence 
class, their inequivalence will most likely be revealed when we calculate their invariants. 

Outline of this article 

We have investigated both extremal and nonextremal PPT states in the system of three qubits using 
both numerical and analytical methods. We classify the states according to the ranks mo, mi, ra-i-, ms 
of p, p Tl , p T2 , and p Ts . 

In Section [2] we present the two main numerical algorithms that we have used in random searches 
for PPT states. We have searched systematically for extremal states of unspecified ranks, and also for 
states of specified ranks that are not necessarily extremal. 

In Section [3] we present results from the searches. We find PPT states with a wide variety of rank 
combinations. The rank 1111 states are the pure states. The states of ranks 2222 and 3333 are all 
separable, in agreement with previously known results. An interesting class of extremal PPT states 
are those of rank 4444, they contain genuine tripartite entanglement, since they are separable in any 
bipartite splitting of the three qubit system. We have studied these states in detail and shown how to 
construct them analytically, as reported in Section |4] 

All the other states we find have all ranks equal to five or higher. We find PPT states, extremal 
or nonextremal, of very nearly every rank combination. For an extremal PPT state the square sum of 
ranks can be at most 193. Up to this upper limit we also find extremal PPT states of very nearly every 
rank combination. 

In Section |4] we present our understanding of the rank 4444 extremal PPT states. It turns out 
that the key to understanding them is the fact that they are biseparable in three different ways lUTI . 
There are two distinct classes of such states, we call them simply type I and type II. The most obvious 
distinction is that the quadratic Lorentz invariant is nonzero for states of type I and zero for states of 
type II. 

When we sort these states further into SL ® SL ® SL equivalence classes, we find that every 
equivalence class of type I states contains a density matrix which is real and symmetric under all 
partial transposistions. We define this matrix to be a standard form for all the states in the equivalence 
class. The UPB states are a subclass of the type I states, but they are not sufficiently generic that we 
find any of them in our random searches. 

It is straightforward to construct the most general density matrix which is biseparable, real and 
symmetric under all partial transpositions. The empirical result of our numerical searches is that this 
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construction reproduces all the rank 4444 PPT states of type I. The SL <S> SL (g> SL equivalence classes 
are parametrized by seven continuous real parameters. 

A rank 4444 PPT state of type II can not in general be transformed to a real form. We have studied 
our numerical examples of such states and observed several special properties that they have. Based 
on these observations we present an analytical construction general enough to reproduce all the type II 
states found numerically. The SL ® SL ® SL equivalence classes in this case are parametrized by one 
continuous complex parameter. 

Our work presented here extends previous studies of PPT states in bipartite composite systems 
|[T8l[T9ll20ll2Tll22ll23l . There are obvious similarities and differences, especially between the 3x3 
and the 2 x 2 x 2 systems. 

2 Numerical methods 

Most of the PPT states that we have studied numerically were found by means of two algorithms that 
search for extremal PPT states or for PPT states with a specified combination of ranks. In this section 
we present briefly these two main algorithms. 

2.1 Random search for extremal PPT states 

Extremal states in V\, the convex set of PPT states of unit trace, are found by an iterative process 
where we start with any PPT state, pick a random search direction restricted to the unique face where 
this state is an interior point, and follow this direction to the edge of the face, which is a face of lower 
dimension. The next search direction is restricted to this new face. The procedure is repeated until a 
face of dimension zero, which is an extremal point of V\ , is located. An extremal point of V\ defines 
an extremal ray, a one dimensional face, of V, the cone of unnormalized PPT states. 

A valid search direction from a given PPT state p G V\ is a nonzero matrix a such that p+ea £ V\ 
for every e in some finite interval, t\ < e < £2, with ei < < £2- This means that a is a traceless 
Hermitian matrix satisfying the four equations 

Pia T tPi = a Ti , i = 0,1,2,3 , (9) 

where Pi is the orthogonal projection onto Img p Ti , and where we write p T ° = p. The real vector 
space H of Hermitian N x N matrices has dimension N 2 and is a Hilbert space with the scalar 
product 

(A,B) = (B,A) =Tr(AB) . (10) 

The constraints in equation (|9) can now be written as eigenvalue equations Pjcr = a for linear pro- 
jection operators Pj on H, symmetric with respect to this scalar product, defined as follows, 

P 4 <7 = (Pi(T Ti Pi) Ti , i = 0,l,2,3. (11) 

The four eigenvalue equations may be combined into one single eigenvalue equation 

3 

^P if T = 4a. (12) 

i=0 

This equation always has the trivial solution a = p, but p is not traceless. If p is the only solution it 
is an extremal point of V\. The number of linearly independent traceless solutions to this eigenvalue 
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equation is the dimension of the restricting face of V\. We find a random search direction by finding a 
complete set of eigenvectors solving equation (fT2l ). then choosing o§ as a random linear combination 
of these, and finally obtaining a traceless solution as 

o = a - (Tr a )p . (13) 

The bipartite version of this algorithm is described in more detail in |[T9l . 

2.1.1 Separability test for low rank states 

A modified version of this algorithm can be used as a separability test for low rank states. If a PPT 
state is not extremal it can be written as a convex combination of extremal PPT states of lower rank. A 
set of such states can be found by searching in both directions a and —a every time a search direction 
a is chosen. The original state is separable if all the extremal states found in this process are pure. In 
general it may be possible to write a separable state as a convex combination involving some mixed 
extremal states. However, for low rank states few mixed extremal states of lower rank are available 
and these "false negatives" are therfore less likely to occur. 

2.1.2 Upper limit on the ranks 

We derive an upper limit on the ranks of extremal PPT states by counting the number of constraints on 
a. Since each operator is a projection with eigenvalues and 1, the constraint equation Pjd = a 
places a number of constraints on a equal to the dimension of the kernel of Pj. If p Ti has rank m«, 
then the rank of Pj is mf, and the dimension of its kernel is N 2 — mf. Thus, for a given p the total 
number of constraints on a, apart from the zero trace condition, is 

4JV 2 -^mf. (14) 

i 

These constraints need not be linearly independent, hence this is an upper limit on the number of 
linearly independent constraints. Since p is an extremal PPT state if and only if it is the only Hermitian 
matrix satisfying the constraints, the number of linearly independent constraints when p is extremal 
must be N 2 — 1. The total number of constraints must be at least as large. This implies the following 
limit on the ranks of extremal PPT states and their partial transposes in a three qubit system, 

rnf < 3N 2 + 1 . (15) 

i 

The derivation of the formula involves a sum over the independent partial transposes, the number 
of which is 2 n_1 for an n-partite system. The formula can therefore be generalized to 

J^mf <(2 n ' 1 -1)N 2 + 1 (16) 

i 

for an n-partite system of total dimension N. For large n this inequality is a very mild restriction on 
the ranks. 
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2.2 Search for PPT states of specified low ranks 



The second algorithm we describe take as input the desired ranks mo, mi,m2,ni3 of p and its partial 
transposes. After we find a p with the specified ranks we may check, by the algorithm just described, 
whether or not it is extremal. 

The equations defining a Hermitian matrix p as a PPT state of the given ranks are of the form 
fj, = 0, where p is a list of all the lowest eigenvalues of p and its partial transposes. We include 
the N — rrii lowest eigenvalues of p Ti (recall our definition p T ° = p). When we write p as a linear 
combination 

j 

of matrices Mj forming a basis for the real vector space of Hermitian matrices, the equations are of 
the form pi(x) = 0. There are N 2 real variables Xj, and the number of equations is 

N e = AN-J2™i- ( 18 ) 

i 

2.2.1 A minimization problem 

One way to solve the equations pi(x) = is to minimize the function 

fix) = J>*(x)) 2 (19) 

i 

and obtain a minimum value of zero. We may use a random search algorithm for finding the minimum. 
This approach may be good enough for many purposes, although not very efficient, and it has the 
advantage of being easy to program. 

Note that a minimum point of a function is in general not very precisely determined numerically, 
because the function varies quadratically around its minimum. This problem does not arise here 
because the desired minimal value is zero. 



2.2.2 The conjugate gradient method 

A more refined and efficient iterative approach is to linearize the equations and solve in each iteration 
the approximate, linearized equations using the conjugate gradient method. A similar algorithm for 
the bipartite case is described in BOl . 

Given an approximate solution x we try to find a better solution x + Ax using the linearized 
equation 

p(x + Ax) = p(x) + B(x) Ax = , (20) 

where the matrix elements of B[x) are Bij(x) = dpi(x)/dxj. We multiply this equation by B T , and 
get another equation 

A Ax = b, (21) 

where A = B T B and b = — B T p. The matrix A is real and symmetric as well as positive semidefinite. 
It is likely to be singular, but the last equation may anyway be solved by the conjugate gradient method. 

We calculate the components of B using first order perturbation theory. If for example pi is the 
eigenvalue of p then we use the formula 

dpi d\ k A dp i Us i 

rtrr j riT ■ *■ rim- • ^ J 
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Here ipk is the eigenvector corresponding to the eigenvalue of p. This formula is based on non- 
degenerate perturbation theory, which is strictly speaking not valid in the present case where we try 
to make many eigenvalues simultaneously equal to zero. The formula nevertheless works well in 
practice. We use similar formulas for the eigenvalues of the partial transposes of p. 

3 Results of numerical searches 
3.1 Ranks of PPT states found numerically 

Numerical searches for PPT states of specified ranks were successful for most of the rank combina- 
tions. Note however that our search algorithm might occasionally produce PPT states of lower ranks 
than the ones specified. Note also that for the purpose of classification of states there is full symme- 
try between p and all its seven partial transposes (including the double partial transposes and the total 
transpose). Furthermore, with three identical subsystems there is full symmetry between states related 
by an interchange of the subsystems. Therefore when we list the ranks mo, mi, mi, 7713 we use the 
convention that mo < mi < 7712 < m.3. 

The PPT states found with specified ranks might be separable or entangled, extremal or nonex- 
tremal. We used them as starting points in searches for extremal PPT states. This would either show 
them to be extremal, or return extremal PPT states of lower ranks. 

A list of all confirmed rank combinations and the corresponding square sums of ranks is given 
in Table [T] By equation ([TBI , for an extremal PPT state the square sum of ranks can not be larger 
than 193. We found no PPT states at all with the very asymmetric rank combinations 5568, 5588, or 
5888. Otherwise, the rank combinations 2222, 3333, and 5688 are the only cases where the inequality 
allows extremal PPT states to exist, but we did not find any. 



1111 


4 


5666 


133 


6677 


170 


2222* 


16 


5667 


146 


6678 


185 


3333* 


36 


5668 


161 


6688 


200 


4444 


64 


5677 


159 


6777 


183 


5555 


100 


5678 


174 


6778 


198 


5556 


111 


5688* 


189 


6788 


213 


5557 


124 


5777 


172 


6888 


228 


5558 


139 


5778 


187 


7777 


196 


5566 


122 


5788 


202 


7778 


211 


5567 


135 


6666 


144 


7788 


226 


5577 


148 


6667 


157 


7888 


241 


5578 


163 


6668 


172 


8888 


256 



Table 1: A list of all confirmed rank combinations and the corresponding sums of squared ranks. The 
upper limit on the sum of squared ranks of extremal PPT states, as given by equation (fT5l) . is 193. 
There are three cases, each marked by an asterisk, where this limit allows extremal states to exist but 
only nonextremal states were found. 

A similar study on the ranks of extremal states in low dimensional bipartite systems is presented in 
GUI . The structure of the rank combinations of our three qubit states closely resembles the structure 
in the 3x3, 3x4, 3x5, and 4x4 systems. For the lowest rank combinations only states where all 
partial transposes have the same rank are found, and below some threshold rank, four in our case, 
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all these states are separable. Above the threshold extremal states are found for essentially all rank 
combinations where PPT states exist and where the upper limit given in equation (fT5l) holds. The sole 
exception is the rank combination 5688 where extremal states are allowed by the limit, but we were 
unable to confirm their existence. We consider it likely that extremal states of this rank combination 
do exist, but are difficult to find numerically by our methods. 

3.2 Rank 4444 PPT states 

The rank four extremal PPT states are of special interest because they are the entangled PPT states of 
lowest rank. Thus, any rank four PPT state is either extremal or separable. Another special property 
of the rank four PPT states is that all their partial transposes have the same rank, we find no PPT states 
of rank 4445, for example. 

Nine out of 196 extremal rank four states found in a random numerical search stand out because 
their quadratic Lorentz invariant is zero, which means in practice of order 1CP 14 to 10~ 16 . We find a 
qualitative difference between the states with nonzero quadratic invariant and those with the invariant 
equal to zero. In particular, in our numerical sample we see no continuous transition from invariants 
of order one to those of order 10 -14 . 

We will discuss the rank 4444 entangled PPT states in much more detail in Section [4] 

3.3 PPT states with special symmetries 

An obvious way to construct a PPT state is to construct a positive Hermitian matrix which is symmetric 
under all three partial transpositions Tf, T2, T3. These symmetries together imply symmetry under the 
total transposition T = T1T2T3, which makes the matrix real. The most general Hermitian matrix 
symmetric under all three partial transpositions contains 27 real parameters. See Appendix lAl 

Completely symmetric PPT states of rank eight are easily generated as follows. Generate a random 
complex matrix, and multiply it with its Hermitian conjugate to obtain a Hermitian and positive matrix. 
Add to this matrix all its seven partial transposes. The resulting matrix is completely symmetric and 
has a high probability of being positive. 

Another way to generate a completely symmetric PPT state of any rank m is the following. Start 
with any matrix in the 27 dimensional space of completely symmetric matrices. Then take the square 
sum of its 8 — m lowest eigenvalues, and minimize this by varying the matrix so as to obtain the 
minimum value of zero. We used this method to create states of rank four, five, six, and seven. The 
rank four states were all extremal, while all the higher rank states were nonextremal. 

It is also possible to relax the symmetry requirements and only require symmetry under two partial 
transpositions, say T\ and T2. This raises the number of free parameters in a Hermitian matrix from 
27 to 36 and allows it to be complex. A state p with the symmetries p = p Tl = p T2 is still guaranteed 
to be a PPT state, because p T:i = p T has the same eigenvalues as p. Only about 75% of the rank four 
states generated with these symmetries were extremal. 

3.4 Decomposition 

The separability test described in Subsection 12.1. H is of particular interest for the nonextremal rank 
5555 states, which are the states of lowest rank where the test is nontrivial. Among the states of this 
kind found in random searches for states of specified ranks we found both separable and entangled 
states. Every entangled state was found to be an interior point of a one dimensional face of V\ , this 
face is then a line segment with a pure state at one end and a rank 4444 extremal PPT state at the 
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other end. Every separable state was found to be an interior point of a four dimensional face of V\, as 
expected when it is a convex combination of five pure states. 

We also tested the rank five states that were constructed so as to be symmetric under all partial 
transpositions. The result is that every such state is an interior point of a face of dimension five, and 
decomposes always into two rank four extremal states having the same symmetry. 

We found the last result surprising, and therefore applied the same test to symmetric states with 
higher ranks. The rank six states behave in a very similar manner, decomposing into four symmetric 
rank four states that are always extremal. The rank seven states, on the other hand, decompose mainly 
into rank 6677 states which obviously do not possess the same symmetiy. Every rank eight state 
decomposes into 32 rank 6777 states that also do not possess any symmetries. 

4 Rank four entangled PPT states as biseparable states 

We find empirically that the entangled and extremal rank 4444 PPT states in dimension 2x2x2 fall 
into two qualitatively different classes, which we choose to call types I and II. By definition, a state p 
of type I has a nonzero (always positive) value of the quadratic Lorentz invariant 

f^px flv = ~Tr{p T EpE), (23) 

whereas this invariant vanishes for a state p of type II. See Appendix |Bj As we shall see, another 
characteristic difference is that a state of type I can be transformed by a product transformation to a 
standard form where it is real and symmetric under all partial transpositions, whereas a state of type II 
can not in general be transformed to a real form. Type I includes as a special case the so called UPB 
states constructed from unextendible product bases (UPBs), introduced in |[T5l[T6l . which we will also 
describe here. 

The basic fact leading to an understanding of these states is that if p is a density matrix of rank 
four in dimension 2x2x2, and if p Tl > 0, then p is separable in dimension 2x4. Similar results 
hold of course for the cases p T ' 2 > and p T ' A > 0. This has been proved analytically ifTTl and also 
verified numerically EOl . We will describe in this section how to use the threefold biseparability in 
order to understand our rank 4444 PPT states. 

We first discuss some common properties of the two types of states, and a new method for con- 
structing them numerically, before we turn to the more detailed mathematical understanding of the 
two types separately, including the UPB states. 

4.1 General considerations 

Consider a generic four dimensional subspace U C <D 8 . It contains exactly four product vectors 

ei = Xi®Ui, i = 1,2,3,4, (24) 

of dimension 2x4, which form a basis for U, and which may be calculated by the method described 
in Appendix |E] It contains also exactly four product vectors 

fj = yj ® s vj , j = 1, 2, 3, 4 , (25) 

where yj G C 2 , Vj 6 <D 4 , and where we use the split tensor product defined in Appendix |F| These 
product vectors form another basis for U. And it contains exactly four product vectors 

9k = Wk ® z k , fe = 1, 2, 3, 4 , (26) 
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of dimension 4x2, forming a third basis for U. 

A density matrix p of rank four with Img p = U and p Tl > must be biseparable in dimension 
2x4, and hence of the form 

4 

p = ^ Ai efej with all A; > . (27) 

i=l 

Note that although p is separable as a bipartite state, it is necessarily entangled as a tripartite state, 
except in the special case when the four vectors ej are 2 x 2 x 2 product vectors, which means that 
they are identical to the vectors fj, and also identical to the vectors g^. 

If p T ' 2 > 0, then p must be biseparable by the split tensor product, and we must have 

4 

P = ^^jfjfj witha11 Mi>0. (28) 
If p T ' A > 0, then p must be biseparable in dimension 4x2, and 

4 

P = ^2 Vk 9k9 k witn a11 u k > . (29) 

fe=l 

4.1.1 Constructing PPT states numerically as biseparable states 

We find that with a generic subspace U any two of the three equations (1271 . (l28l i. and (1291 are incom- 
patible, they have no common nonzero solution for p. If we require that two of the three equations, 
for example (l2~7T ) and (l28l i. should be compatible, then this restricts the subspace U. It is easy to find 
subspaces where two of the three equations are compatible, but then the third equation will in general 
be incompatible with the two. This means that we may have, for example, p > 0, p Tl > 0, and 
p T ' 2 > 0, but p T:i ^ 0. 

When we want a PPT state with p > and p T% > for i = 1, 2, 3, we need to find a subspace U 
where the equations (1271 and (|28T ) are compatible, and at the same time the equations (T27T ) and (1291 . 
We will describe next how to solve these compatibility problems numerically. 

The equations (1271 and (|28T ). apart from the positivity conditions on the coefficients Aj and pj, are 
compatible when there is at least one linear dependence between the eight Hermitian 8x8 matrices 
e^e] and fjfj- We do the calculation numerically by converting the eight matrices to eight vectors in 
IR 64 and making a singular value decomposition of the 64 x 8 matrix. The eight singular values are 
non-negative, by definition, and we want the smallest one to be zero. This is a minimization problem 
where we vary the subspace U in order to minimize the smallest singular value. We solve this in 
practice by a random search method. 

An output of the singular value decomposition, after minimization, is the two sets of coefficients 
Xi and pj corresponding to the singular value zero. However, the singular value decomposition puts 
no restriction on the signs of the coefficients A» and pj. If we get a solution with all Aj < 0, then we 
simply switch all four signs in order to get all Aj > 0. If, on the other hand, we get four coefficients 
Xi with both signs, then we have a matrix p which is biseparable in two ways, but we can not make it 
positive semidefinite. This is then a subspace U which we can not use. 

The procedure for making the equations (l2~7T i and (|29l compatible is the same, we formulate it as 
a second minimization problem for a smallest singular value. In order to obtain a subspace U solving 
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both problems simultaneously we simply minimize the sum of the smallest singular value from the 
first problem and the smallest singular value from the second problem. 

The numerical procedure described here for constructing PPT states of rank four works well, 
although the random search for a minimum may take some time. In practice, we always get one 
singular value zero with a corresponding unique solution for the coefficients Aj, and z/&. However, 
it seems that each coefficient Aj for % = 1, 2, 3, 4 comes with an essentially random sign. We need the 
same sign for all four coefficients in order to satisfy the positivity condition p > 0, and on the average 
we have to try about eight times before we succeed once. 

The majority of the rank four PPT states found numerically by this method are of type I, having 
a nonzero value of the quadratic invariant. We find however also a significant fraction of states of 
type II, with vanishing quadratic invariant. As already remarked, these are two distinct classes of rank 
four PPT states. We will proceed next to discuss them separately, after discussing the so called UPB 
states. 

4.1.2 Standard form 

In the generic case there always exists a product transformation 



uniquely defined up to normalization, transforming the two dimensional vectors x, y, z into V\x, V2V, 
V3Z such that the transformed vectors after suitable normalizations have the standard form defined in 
Appendix O equation (1 100b . with complex parameters ti, t2,t$, 



If p is a density matrix of rank four with Img p = U, then the product transformation V transforms p 
into VpVK which in our terminology is SL (g) SL <g> SL equivalent to p. This transformed matrix will 
then be a standard form for p. 

Note that the standard form is not completely unique, because it depends on the ordering within 
each of the three sets of product vectors ei, fj, and g^. Thus there is a discrete ambiguity. 

4.2 States constructed from unextendible product bases 

The construction of rank four entangled PPT states in dimension 2x2x2 from unextendible product 
bases (UPBs) is due to Bennett et al. fl~5l [Toll . It has been shown, both numerically ETl 1231 and 
analytically |[24ll25l . that the UPB construction in dimension 3x3 produces all the rank four entangled 
PPT states there. In contrast, we find that the UPB states in dimension 2 x 2 x 2 are only a subclass 
of the rank four PPT states of type I, with nonvanishing quadratic invariant. 

An unextendible product basis (a UPB) is a maximal set of orthogonal product vectors which is 
not a complete basis for the vector space. In other words, it is an orthogonal basis of product vectors 
for a subspace, with the property that the orthogonal subspace contains no product vectors. As shown 
in Appendix iDl a UPB in dimension 2x2x2 consists of four product vectors 



V = Vi®V2®V 3 



(30) 




(31) 



ipi = Xi ® yi ® Zi 



i = 1,2,3,4 . 



(32) 
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From these vectors, normalized to unit length, we construct a PPT state p as a normalized projection 
operator onto the orthogonal subspace, 

^ i=l ' 

where 1 is the 8x8 unit matrix. This is obviously a PPT state, since p Tl , for example, is again an 
orthogonal projection, 

^ i=i ' 

with 

i>% = x* ®yi®Zi. (35) 

The state p constructed in this way is obviously entangled, since Img p contains no product vector, 
and a characteristic property of a separable state a is that Img a is spanned by product vectors. From 
p we get other rank four entangled PPT states of the form p = aVpV^ where V is a product matrix, 
V = V\ ® V2 <8> V3, and a is a normalization constant. We have then that 

Imgp = VImgp, Kerp = (V^Kerp . (36) 

The kernel of p is spanned by the four orthogonal product vectors fy. The kernel of p is spanned by 
the product vectors 

(^rv* = (o^r 1 ^) ® (o^r 1 ^) ® (c^r 1 ^) > 

which are orthogonal if and only if V is unitary. 

The existence of a basis of product vectors in the kernel, orthogonal or not, is the characteristic 
property of the rank four PPT states of the UPB type and their SL ® SL (g> SL transforms. These states 
are not generic, because a generic subspace of dimension four contains no product vector. In fact, a 
generic and suitably normalized product vector 




contains three complex parameters a, b, c. In order to lie in a given four dimensional subspace it has 
to satisfy 8 — 4 = 4 linear constraints, which it can not do in general with only three parameters. 

Using numerical methods, we have searched for product vectors in the kernels of our extremal 
rank four PPT states, but found none. We conclude that there are no states equivalent to UPB states 
in our numerical sample. The UPB states exist, but are not sufficiently generic to be easily found in 
random numerical searches. 

The UPB has a standard form where all the vectors x, y, z are real and given by three continuous 
real parameters (three angles), as in equation (11091 . Therefore p has a standard form in which it is 
symmetric under all partial transpositions, p = p Tl = p T-2 = p T:i . 

4.3 Rank four PPT states with nonzero quadratic invariant 

We turn now from the UPB states to the more general class of rank 4444 extremal PPT states with 
nonvanishing quadratic invariant. The key to understanding these states comes from a study of our 
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numerical examples. It turns out that all the states we have found numerically may be transformed to 
a standard form where they are symmetric under all partial transformations, and hence real. 

Let p be a state with nonvanishing quadratic invariant. We compute the three sets of four product 
vectors ej = x\ <g> U{, fj = yj® s Vj, and gk = Wk® Zk in Img p, with Xj, yj,Zk £ <D 2 and Ui,Vj,Wk G 
<D , as already defined in Subsection 14.11 Then we use a product transformation V = V\ ® V\ <g) V3 to 
transform the two dimensional vectors x, y, z to the standard form given in equation (1311 . 

We find empirically that we always get real values for the parameters ti, t2, £3, and that the same 
transformation V always makes u, v, w real. Moreover, the transformed matrix VpV^ always becomes 
symmetric under all partial transpositions. These are the key observations. We have no mathematical 
proof that p must have these properties if its quadratic invariant is nonzero, but we take it as a working 
hypothesis. As such it is very powerful, and we will see next how to construct the most general PPT 
state with these properties. 

4.3.1 Explicit construction 

In equation (1271) we may absorb a factor y/Xi into the vector = X{ <X> U{ and write 

4 

p = 5>et. (39) 

i=l 

The problem we face is to construct a density matrix p which has this form with x and u real, and is 
symmetric under all partial transpositions, p = p Tl = p T2 = p Ti . We also want x to have the standard 
form given in equation (DTI ), so that 

e 

Then p has the following form, 



with 

A 
B 

C 

Let us see what happens if we simply take u to be a random 4x4 real matrix. As shown in 
Appendix |Aj p as given in equation (I4TI) is symmetric under all partial transpositions if and only if the 
4x4 matrix A is real and has the following general form, 

/ ai a 5 a 6 
j, ^ a 5 a 2 a 7 
a 6 a 7 a 3 
\a 7 a 8 a g 

and B and C are also real and have similar forms. In the case of A the only condition which is not 
automatically satisfied as a consequence of equation (l42l is that An = For B and C there are 
similar conditions. In the case of B the condition takes the form 

- W43U13 + tiUuUU = -U33U23 + il«34«24 ■ (44) 



m 

u 2 



U 3 tlU4 

-u 3 n 4 



(40) 



A B 
B C 



(41) 



U2U2 + U3U3 + u 4 uj . 



(42) 



a 7 \ 

a 8 
a 9 

04/ 



(43) 
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(45) 



(46) 



We solve it simply by choosing 



v%(e <g> e) u 3 

H — —T-, ; — • 

u\ (e <8 e) U4 

Next, the condition that ensures the right form for A is the following, 



(47) 



Ui(e <8 e) u\ + uj(e <8 e) u% + uj(e <8 e) u± = . 



(48) 



Define now 



(49) 



If a 2 > 0, then we replace u\ by mil and have a solution of equation (1481 . If instead a 2 < 0, then 
we change the sign of uf(e (8) e) u\, for example by changing the signs of the first two components of 
u\. With this new vector u\ we have a 2 > 0, and again we solve equation (|48T ) by substituting au\ 
for tii. 

The condition that ensures the right form for C is the following, 



We solve it for U2 in the same way as we solved equation (1481 for u\. This completes the construction 
of p, except that we should normalize p to have unit trace. 

4.3.2 Parameter counting 

A random 4x4 real matrix u contains 16 freely variable real parameters. In the process of constructing 
p to be symmetric under partial transformations we had to redefine the normalizations of u\ and U2, 
this reduces the number of free parameters to 14. The final normalization of p reduces the number of 
parameters to 13. 

Note that when we take x to have its standard form and choose u randomly, the resulting vectors y 
and z will most likely not have their standard forms. In order to transform them to standard forms we 
will need two SL(2, R) transformations, containing altogether six real parameters. After we fix the 
standard forms of x, y, z there are no continuous degrees of freedom left in the product transformation 
V. There is however some discrete freedom left, because we may permute the product vectors ej, as 
well as the vectors fj and the vectors g^. 

We conclude that there is a family of SL (8 SL % SL equivalence classes of rank 4444 extremal 
PPT states with nonzero quadratic invariant described by 13 — 6 = 7 continuous real parameters. 

4.4 States with vanishing quadratic invariant 

In our numerical random searches for PPT states of rank four we found a special class of rank 4444 
PPT states with the property that the quadratic Lorentz invariant of each state p vanishes, that is, 



«^(e (8 e) u 2 + (e <8> e) u 3 + uj (e (8 e) U4 = . 



(50) 



Tr(p J EpE) = . 



(51) 
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See Appendix iBl In this subsection we present our understanding of these states. We find that such a 
state has a standard form described by one continuous complex parameter. 
The spectral representation of p, 

4 

1=1 

with four eigenvalues Aj > and corresponding orthonormal eigenvectors rji, gives that 
Tr(p T EpE) = Y^K^TrlvtvfEvjrijE] = £ A.A, [vl Erjj] [v}Er,*] 

= -J2 XiXj [rjfEty] [4E^} =-J2 X t X, \vfErjjl 2 . (53) 

Since the four eigenvalues are positive, the only way to have Tr(p T EpE) = is to have all the 
eigenvectors orthogonal in the scalar product defined by E, 

vfErt^O for i,j = 1,2,3,4. (54) 

Since E T = — E, this scalar product is antisymmetric, ip T E(fi = —4> T Eip, hence every vector is 
orthogonal to itself, ip T Eip = 0. 

Equation (1541 ) means that the antisymmetric scalar product vanishes identically in the four dimen- 
sional subspace Img p spanned by the eigenvectors rji with i = 1, 2, 3, 4. We may use this observation 
in order to search numerically for PPT states of this type. 



4.4.1 A random search method 

It is a straightforward exercise to construct a random four dimensional subspace U C <D 8 where the 
antisymmetric scalar product vanishes identically. We may start with a random normalized vector rp\. 

t T 

The two conditions ip[ip = 0, ipf Eip = restrict tp to a six dimensional subspace, and we choose 
^2 as a random normalized vector in this subspace. Next, the conditions ipjip = 0, ipfEip = for 
i = 1, 2 restrict tp to a four dimensional subspace, and we choose tp^ as a random normalized vector 
in this subspace. Finally, we choose V>4 as a random normalized vector in a two dimensional subspace. 
In this way we get vectors ifrx , V>2 , ^3 , 4>4 such that 

= Sij , iifEiPj = . (55) 

The vectors <pi = Eip* for i = 1, 2, 3, 4 lie in the orthogonal subspace U^, since 

<t>% = -VfEtyj = , (56) 
and they satisfy the same relations as the vectors ipi, 

4<t>j = -tfE 2 r 3 = (4^)* = % , 

tfEfa = -4E^* = tyjE^y = . (57) 

For a given subspace U we may try to construct numerically a PPT state p with Img p = U,by 
the method described in Subsection l4.1.1l Again we find that the construction fails in general, but we 
may select a suitable subspace with vanishing antisymmetric scalar product where the construction 
succeeds. 
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As remarked in Subsection 14.1.11 when we search for a random subspace U where we can find a 
matrix p satisfying all three biseparability conditions (1271 ). (|28T ). and (1291 ), disregarding the positivity 
conditions, then it happens only in about one case out of eight that the p we find is positive. Sur- 
prisingly, when we now do the same search restricted to subspaces with the special property that the 
antisymmetric scalar product vanishes identically, then the positivity conditions hold, not every time, 
but almost every time. 

4.4.2 Explicit construction 

By studying our numerical examples of such states we find empirically that they have certain proper- 
ties which enable us to construct them explicitly. Again we have no strict proof that these properties 
are necessary, but it is a very powerful working hypothesis to assume that they hold. 

Like before, we introduce three sets of basis vectors for the subspace Img p, ei = xi ® ui, fj = 
yj <g> s vj, and g k = w k ®z k , with x i: yj,z k G <D 2 and Ui,Vj,w k G <D 4 . 

The vanishing of the antisymmetric scalar product in Img p means that 

efEej = (xjex j )(uj(e®e)u j ) = for i,j = 1,2,3,4 . (58) 
We have that xfexj = for i = j, but xfexj / for i ^ j. Hence we must have 

uj (e ® e) Uj = for i^j. (59) 
Note that e is antisymmetric, but e ® e is symmetric, 

(e<g>e) T = e T Oe T = (-e) <g> (-e) = e . (60) 

In our numerical examples we see that equation (1591 is solved in the following remarkable way. 
The vectors Uj with i = 1,2,3,4 can not all be product vectors, because then p would be separable. 
Instead, each m is a linear combination of two product vectors, as follows, 

Ui = a ik i y k ® zi + a imn y m ® z n . (61) 

There is no sum here over the indices k,l,m,n, and we have always k ^ m and I ^ n. For each 
value of i, six different combinations occur for klmn, and we may order the vectors y and z in such a 
way that we get the index combinations listed in Table |2] Since 

(Vk 8> zi) T {e ® e) (y m ® z n ) = (y{ ey m ) (zfez n ) = if k = m or l = n , (62) 

we are guaranteed that m and Uj satisfy the orthogonality relation in equation (l59l if 

Ui = a ik i y k ® zi + a imn y m ® z n , 

uj = a jkn y k ®z n + dimi y m ® z\ . (63) 

We see from Table |2]that two vectors u\ and Uj are always orthogonal in two ways. For example, 

u\ = a lw yi®z 2 + a 12 i y 2 ®zi = a lu y 3 ® z 4 + ai 43 y 4 <8> Z3 , 

u 2 = a 2 n yi ® z 1 + a 22 2 2/2 <8> z 2 = a 23 3 2/3 ® z 3 + a 2 44 2/4 ® Zi . (64) 

We now transform x, y, z to the standard form defined in equation (BTT l. with complex parameters 
t±,t2,t 3 . We find that the linear dependencies listed in Table |2] require that t% = £3, and when this 
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i 






kl 


ran 






1 


1221 


1331 


1441 


2332 


2442 


3443 


2 


1122 


1342 


1432 


2341 


2431 


3344 


3 


1133 


1243 


1423 


2134 


2244 


3241 


4 


1144 


1234 


1324 


2143 


2233 


3142 



Table 2: The allowed index combinations in equation (I6TT) . 



relation holds they give a unique solution for u depending on the single complex parameter t = ti 
t%. The solution is 

f t t t \ 
1 1 -t 
-10 1 -t 

V o i-i -V 

The condition that the four vectors = Wf. ® must be linear combinations of the vectors e 
Xi <X> Ui requires that also t\ = t. We find the overall solution v = w = u, giving the vectors 



u 



(65) 
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* 2 \ 
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1 
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-t 
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-t 
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-i 2 
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-1 


-t 










-1 


-t 
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-t 
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-1 
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-t 2 
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-t 2 








-1 


-t 










-1 


-t 










-1 


-t 








-1 


-t 










-1 


-t 










-1 


-t 


\o 


1 


1 


-1/ 




\o 


1 


1 


-1/ 




Vo 


1 


1 


-1/ 



(66) 



Thus we arrive at the following explicit standard form for p, depending on the single complex 
parameter t, 



9i9, 



i=l 



with 



and 



\t\ 2 \l + t\ 2 



i=i 



i=i 



5|i| 4 + 10|t| 2 + 1 + (3|t| 2 + 1) |1 + 1\ 2 



(67) 

(68) 
(69) 



4.4.3 The effect of partial transposition 

The effect of the partial transposition T\, for example, on the PPT state p given in equation (1671 ) is that 

l 



P T1 



i=i 



(70) 



with e"j = x* (8) ttj. The two dimensional vectors x are complex conjugated while the four dimensional 
vectors n are unchanged. If x and u have their standard forms as given in equation (l3~0 (with ti = t) 
and in equation (I65T ). it means that the parameter f is complex conjugated in x but not in u. 
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We find that the complex conjugation of the standard form of u is equivalent to a linear transfor- 
mation on product form. That is, there exists a 2 x 2 matrix W such that 

{W ®W)u i = C i u* , (71) 

with four complex normalization constants Q having different phases but equal absolute values, 
\Ci\ = \Cj\ for i,j = 1, 2, 3, 4. There are four solutions for W, since the general solution 

w_ f-eit*(l-e 1 \t\+e 2 \l + t\) \t\ (t* + e x \t\) \ 

V |*| + ei** \t\(l-e 1 \t\ + e 2 \l + t\)J { ' L) 

contains two arbitrary signs ei = ±1 and e 2 = ±1. 

It follows when we define V = 1 <g> W W, with the 2 x 2 unit matrix 1, and introduce a 
normalization factor b, that 

bV P Tl V ] = p* = P T . (73) 

This shows that the partial transpose p Tl has a standard form like the standard form of p, but with the 
complex conjugated parameter value t*. The same observation applies of course also to the partial 
transposes p T ' 2 and p Ts , and to the total transpose p T = p* . 



5 Summary and outlook 

We have presented here a numerical survey of entangled PPT states in the system of three qubits, with 
an emphasis on the extremal PPT states. Important characteristics of a state are the ranks of the state 
itself and of its three single partial transposes. 

For ranks equal to five or higher very few rank combinations are missing, the variety could not be 
much larger than what we see. This is true also for the extremal PPT states, where there is an upper 
bound of 193 for the square sum of ranks. The existence of a large variety of very different states 
indicates that it will not be easy to understand all the extremal PPT states analytically, in the same 
way as we understand the pure states. 

A natural place to start a project to gain some analytical understanding is with the rank 4444 states, 
which are the lowest rank entangled PPT states. Another good reason for our interest in these states is 
that they have genuine tripartite entanglement, since they are separable with respect to any bipartition. 
An unanswered question about the rank four states is why the rank combination 4444 is the only one 
observed. 

The most important advance reported here is that we have uncovered the mathematical structure 
behind the rank 4444 entangled PPT states, and shown how to reproduce analytically all the states 
found numerically. They fall into two distinct classes with very different analytical representations. 
A sobering fact is that the known states of the UPB type were not found numerically, although they 
appear as just a special case in our classification scheme. We take this as a reminder that there may 
exist unknown types of rank four PPT states that we have missed because they are not generic. 

We intend to continue our project and try to understand for example the rank 5555 extremal PPT 
states in the three qubit system. The numerical study of higher dimensional tripartite systems quickly 
becomes impractical, simply because the dimensions grow rapidly. 
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A Partial transpositions 



Write an 8 x 8 complex matrix as a 4 x 4 matrix of 2 x 2 matrices, 



X = 



( A B C D\ 

E F G H 

I J K L 

\M NO Pj 



The partial transposition T\ moves 4x4 submatrices, 



(A B I J\ 

E F M N 

C D K L 

\G H O PJ 



T2 moves 2x2 submatrices within 4x4 submatrices, 



X T2 = 



(A E C G\ 

B F D H 

I M K O 

\j N L Pj 



T3 transposes the 2 x 2 submatrices, 



X T 3 



(A T B T C T D T \ 

E T F T G T H T 

I T J T K T L T 

\M T N T O t P T 



J 



(74) 



(75) 



(76) 



(77) 



If X is symmetric under all three partial transpositions, then it has the form 



X 



/A 
B 
C 

\D 



B 
F 
D 
H 



C 
D 
K 
L 



H 
L 



with A T = A, B T = B, and so on. In particular, X is symmetric, X 
addition X is Hermitian, then it is real. 



(78) 



B SL(2, C), Lorentz transformations, and Lorentz invariants 

Let e be the two dimensional Levi-Civita symbol, 

<-(-°. 

It is a square root of —1, e 2 = —1. The Lie group SL(2, C) consists of all 2 x 2 complex matrices 
with unit determinant. If V € SL(2, C), 
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then 



V = 



d -b 
-c a 



-eV 1 e . 



(81) 



Thus VeV T = e, in this sense e is an invariant tensor under SL(2, <D) transformations. 
A general 2x2 Hermitian matrix may be written as 



X = 



n o i 

X + X x - 

1 -9 

X + IX X 



IX 



arcr 



n > 



(82) 



where x^ is a real fourvector, <7o = 1 is the unit matrix, and Oj for j = 1, 2, 3 are the Pauli matrices. 
The determinant of X is 



det(X) 



1 



Tr(X T eXe) = g^afx" , 



(83) 



where is the metric tensor, 



(\ 
0-10 

0-1 

\0 






-v 



(84) 



More generally, for X = x^a^ and Y = y^a^ we have 

Tr(X T eYe) = -2^„zV = -2^^ . 



(85) 



The transformation X >->■ X = VXV^ with det V = 1 is a continuous Lorentz transformation. 
It leaves the determinant invariant, and leaves the scalar product between two fourvectors invariant 
because det V T = det V* = 1, hence V T eV = VUV* = e and 



Tx(X T eYe) = Tx{X T {V T eV)Y (V ] eV*)) = Ti(X T eYe) . 



(86) 



The parity inversion x 2 = —x 2 takes the form X = X T and leaves the scalar product invariant, 
although it is not of the form X = VXV^. 

In <D 8 = <D 2 <g> C 2 (g> C 2 the antisymmetric tensor 



E = e< 



e = 























1\ 




















-1 




















-1 




















1 




















-1 




















1 




















1 




















V-i 




















0/ 



(87) 



is invariant under SLg)SL(X>SL transformations, in the sense that VEV T = E when V = V\®V2®V^ 

and V^V&.VSj G SL(2,C). 

A general 8x8 Hermitian matrix may be written as 



(88) 
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with 4 x 4 x 4 = 64 real coefficients 

a i»>\ = Xr(A ( (g, av ® CTa) ) . (89) 

o 

A product transformation of the form A = VAV^ with V = V\ <8> V2 <8> V3, as above, acts as three 
independent continuous Lorentz transformations on the three Lorentz indices. Note that the partial 
transpositions are discrete Lorentz transformations, since they are parity inversions. 
For Hermitian matrices A, B the quantity 

Tt(A t EBE) = -8g Xa g^g^a x ^b a ^ = -8a x ^b x ^ (90) 

is real and invariant under the product transformations A = VAV^, B = VBV^ . It is also invariant 
under all three partial transpositions. Note that each Lorentz index on a tensor a x ^ u represents its own 
subsystem and can therefore only be contracted against the corresponding index on the tensor b x ^ v '. 
A density matrix p in dimension 2x2x2 has one quadratic Lorentz invariant 

h = P^P^x = ~\ Tr(p T EpE) > . (91) 
For pure states pi = ipiipj we have that 

Tr(pjE Pj E) = [iijE^li^Etf] = -[^E^][4 E ^] = -|^^f • (92) 
The inequality I2 > follows because p is always a convex combination of pure states, 

p = Y J \^i4 with A. t >0, X)Ai = l, (93) 

i i 

and hence 

Tr(p T EpE) = - Ai Xj I ifiEfy 1 2 < . (94) 

i,j 

There are five different fourth order invariants obtained by different combinations of index con- 
tractions, but one of these is simply the square of the second order invariant. The four new invariants 
can be written as 

hi = P^P^P^PaVX, 

142 = ^ Ptf\ P^ 1 Pau-y , 

143 = P 11 ^ P^ P^ PauX , (95) 
I U = P^p/ 7 P\yPapX. 

Note that all these Lorentz invariants are invariant under SL(g)SL(g>SL transformations of a density 
matrix without subsequent normalization to unit trace. Division of the fourth order invariants by the 
square of the second order invariant gives true invariants that are also independent of the normalization 
of the density matrix. They may be used in order to test whether two density matrices belong to the 
same SL (g> SL ® SL equivalence class. 
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C Standard forms of sets of vectors 

Given four vectors X{ G <D 2 , i = 1,2, 3, 4. We write 

x = (x 1 ,x 2 ,x 3 ,x 4 ) = C d 6 f 9 X (96) 
We consider here the generic case with det(x{, Xj) / for % ^ j. Multiplication by the matrix 

and a subsequent normalization of the vectors gives the form 

y= Vo 1 h i)> (98) 

with det(yi,yj) ^ for j 7^ j, which means that £1*2 7^ 0, 1. Multiplication by 

V-(-* J) (99) 



and normalization now gives the standard form 



1 1 t 
1 -1 1/' 



(100) 



with one variable parameter t = — M2 7^ 0, —1. Since 



t _ det(zi,z 3 )det(z 2 ,^4) 
det(^i, za) det(z 2 , -23) 

and this ratio of determinants is invariant under nonsingular linear transformations and vector normal- 
izations, we have that 

t= det(si, S3 )det (3:2,3:4) = ( a f ~ be)(ch - dg) 
det(xi, X4) det(x2, 3:3) {ah — bg)(cf — de) 

We see that t will be complex in the generic case. In the special case where t is real and positive, we 
may multiply by 

W= I V\ (103) 



v o Vt, 

and normalize so as to obtain the standard form 

w= (o 1 -St ij' (104) 

where the vectors are real and pairwise orthogonal, wjwj = wfwj = for i, j = 1, 2 and i, j = 3, 4. 
Instead of equation (1 1001) we might have chosen one of two alternative standard forms, 
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or 



1 1 if' 
0-111 



The invariant formula for t, equation (1 100b . gives that 



or inversely, 



t' 



■1-t', 



1 - t 



1 



1 + t" 



(106) 



(107) 



(108) 



In the case where t is real we see that t > gives t' < 0, t" < 0, whereas — 1 < t < gives t' < 0, 
t" > 0, and t < — 1 gives t' > 0, t" < 0. Thus, if t is real and t ^ 0,-1, there is always exactly 
one pairing of the four vectors, either (X1X2) (X3X4), (X1X3XX2X4), or (xiX4)(x 2 X3), such that there 
exists a linear transformation which will make both pairs real and orthogonal. 



D Unextendible product bases 

An unextendible product basis (a UPB) in a subspace U of <D 2 C 2 £8) <D 2 is a set of orthogonal 
product vectors e« = Xi yi ® z% spanning U such that no product vector is orthogonal to them all. 
We assume that the vectors Xi, Xj £ <C 2 are linearly independent when i 7^ j, and similarly with yi, yj 

and Zi , Zj . 

Obviously, in order to have e, _L ej we must have either Xi _L Xj, yi _L yj, or Zj _L zj. These 
conditions have a solution with four product vectors, unique up to permutations, as illustrated in 
Figured] In a similar way as discussed in Appendix ICl the vectors may always be transformed to the 
real standard forms 



x 



1 
1 



Si Cl 



c 2 s 2 
-s 2 1 c 2 



C3 s 3 

-S3 C 3 1 



(109) 



where c, = cos Sj = sin 6i, and 61,62, 63 are three angular parameters. 




Figure 1: Orthogonality relations of an unextendible product basis = Xj (g> yi ® Zj, i = 1, 2, 3, 4. 
An unmarked line: Xj _L Xj. A line with one tick mark: yi _L yj. A line with two tick marks: Zj _L Zj. 
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E A problem of finding product vectors in a subspace 

1, . . . , 8, j = 1, . . . , 4. We write a linear 



a 



Given four vectors ipj 6 <D with components tpij, i 
combination of them as a matrix product, 

4 
3=1 

with a € <D 4 . Assume that <fi is a tensor product 



/ oce\ 
ac/ 
a<ie 

;K)*C)- - 

6c/ 
bde 
\bdfj 

The presence of the first factor in the tensor product implies the equation 

(A - [xB)a = , 
where fi = a/b, and A and B are the following 4x4 matrices, 





flpll 012 


013 


0m\ 








052 


053 


^54\ 


A = 


021 ^22 


023 


V>24 


B 




061 


^62 


^63 


^64 


031 ^32 


^33 


V>34 




071 


^72 


^73 


^74 




\1p41 ^42 


^43 


044/ 






W81 


^82 


^83 


084/ 


The presence of the second factor implies the equation 


















(C- 


- /iD)a = 


o, 










where /j, = c/d, and where 




















^011 012 


013 








/031 


032 


^33 


^34 \ 


C = 


021 022 


023 




D 




041 


V>42 


^43 


^44 


051 052 


^53 


V>54 




071 


V>72 


073 


t/j U 




V061 ^62 


^63 


^64/ 






V081 


^82 


^83 


1p84/ 


The presence of the third factor implies the equation 


















(E- 


- fiF)a = 


o, 










where fi = e/f, and where 




















/011 012 


013 








/0>21 


022 


023 


024\ 


E = 


031 032 


^33 


^34 


F 




041 


^42 


V^43 


^44 


051 052 


053 


^54 






061 


^62 


^63 


^64 




V071 V>72 


^73 








\081 


^82 


^83 


084/ 



(110) 



(111) 



(112) 



(113) 



(114) 



(115) 



(116) 



(117) 
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Equation (II 121 ) is a generalized eigenvalue equation, having in the generic case four different 
complex eigenvalues \ii with corresponding eigenvectors on, defining four vectors that are product 
vectors in dimension 2 x 4 of the form 



(pi = ipai = Xi ® Ui with Xi= yiJ> Ui = Bai. (118) 

Similarly, equation (11141 ) gives four eigenvalues ^i with corresponding eigenvectors a-i, denning 
four vectors that are product vectors when we use the split tensor product denned in Appendix 10 

(pi = ipai = yi® s Vi with yi = yf) > Vi = Doii ' ^ 119 ^ 

Finally, equation (11161 ) gives four eigenvalues )Uj with corresponding eigenvectors ai, denning 
four vectors that are product vectors in dimension 4x2, 

(pi = ipai = Wi® Zi with Z{ = , Wi = Fai ■ (120) 

If the two equations (11121 ) and (11161) have a common eigenvector a, then <p = tpa is a product 
vector in two ways, both in dimension 2x4 and in 4 x 2. This means that it is a product vector in 
dimension 2x2x2, and the same a is an eigenvector of equation (II 141 ). 



Note that the standard method for solving the generalized eigenvalue equations (11121) . (1 1 14b . 
and (11161 ) depends on the nonsingularity of the matrices B, D, and F. If one or more of these matrices 
are singular, we may usually avoid the problem simply by making a random product transformation 
ip i — y ip = Vip with V = V\ ® V2 8) V3, then solving the problem with tp instead of ip and transforming 
back in the end. 



F The split tensor product 



We find it useful to define a split tensor product so as to be able to take out the middle factor in a 
tensor product of three factors. Thus we define 



x ® y ® z = y ® s (x <2> z) . 
For y = (c, d) T <E <D 2 and v = (p, q, r, s) T G C 4 = C 2 ® C 2 we define 



fcp\ 
cq 
dp 
dq 
or- 
es 
dr 

\dsj 



y® s v 



(121) 



(122) 



This corresponds to equation (lllll i with p = ae, q = af, r = be, s = bf. 
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